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Abstract 

We consider the Cauchy problem for the Perona-Malik equation 

Uf = div ,„ lo 

VI + |Vw| 2 

in an open set Q C R n , with Neumann boundary conditions. 

It is well known that in the one-dimensional case this problem does not admit any 
global C 1 solution if the initial condition uq is transcritical, namely when |Vmo(#)| — 1 
is a sign changing function in Q. 

In this paper we show that this result cannot be extended to higher dimension. We 
show indeed that for n > 2 the problem admits radial solutions of class C 2,1 with a 
transcritical initial condition. 
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1 Introduction 



Let O C 1" be an open set. Let us consider the Cauchy boundary value problem 

V(x,t) G n x [0,T), (1.1) 

V(x, t) G dQ x [0, T), (1.2) 
«o(a;) G a (1.3) 

This problem was considered in [17] in the context of image denoising. In that 
framework uq is the grey level of a (noisy) picture defined in a rectangle Q C IR 2 , and 
u(x,t) for a given t > should represent a restored image obtained by smoothing the 
regions where |Vw | < 1 (objects) and enhancing the regions where |Vw | > 1 (edges). 

Equation (1.1) is the formal gradient flow of the functional 

PM{u) := \ I log (1 + \Vu{x)\ 2 ) dx. 
2 Jn 

The convex-concave behavior of the integrand makes (1.1) a forward-backward par- 
tial differential equation of parabolic type, with a forward (or subcritical) region where 
\Vu(x, t)\ < 1, and a backward (or supercritical) region where \Vu(x, t)\ > 1. If the 
initial condition is subcritical, namely |V-u (a;)| < 1 for every x G Q, then the maximum 
principle guarantees that the same is true for all subsequent times (see [14]). It fol- 
lows that in this case equation (1.1) is always forward parabolic, and a smooth solution 
globally exists (T = +oo). 

Things are far more complicated when the initial condition is trans critical, namely 
when there are points x G Q where |Vuo(^)| < 1, and points x G f2 where |V«o(^)| > 
1. In this case the forward-backward character of the equation makes the problem 
ill posed from the analytic point of view (see [15]). On the other hand, numerical 
computations exhibit much more stability and show that the process has the desired 
denoising effect on the initial condition. This is usually referred as the Perona-Malik 
paradox. The mathematical understanding of this phenomenology is still far away, 
despite of the considerable amount of interest generated by this problem in the last 
fifteen years (see [2, 3, 4, 5, 6, 7, 8, 12, 13, 18, 19]). 

In this paper we focus on classical solutions, namely solutions which are at least 
of class C 1 or C 2 ' 1 . Throughout this paper C 2,1 denotes the standard parabolic space 
of functions with one continuous derivative with respect to time, and two continuous 
derivatives with respect to space variables. 

As far as we know, in literature there are two main results for classical solutions, and 
they are both limited to the one-dimensional case. The first known result was proved 
in [10] and concerns local-in-time solutions. The result is that for n — 1 the set of 



/ \ ,. / Vu(x,t) \ 

du , 

Tn^ = 
u(x, 0) = 



1 



initial data for which a local-in-time C 2,1 solution exists is dense in C 1 (r2). On the other 
hand, one cannot expect existence of local classical solutions but for a dense set of initial 
conditions (see [15]). 

The second known result concerns global-in-time solutions, and it is a nonexistence 
result. It states that when n = 1 there does not exist any global-in-time solution of 
class C 1 with a transcritical initial condition. This was first proved in [14] with some 
technical assumptions on u , and then in [11] in full generality. 

Many attempts have been made to extend this result in higher dimension. The main 
point in the one-dimensional proof is the so called persistence of supercritical regions, 
namely the fact that supercritical regions, if present for t — 0, cannot disappear for all 
subsequent times. This important qualitative property of solutions is also consistent 
with intuition, because it means that the process does not destroy edges. 

In dimension n — 1 the persistence of supercritical regions follows from the fact that 
the L°° norm of |Vw| in Q is an increasing function of time. Unfortunately this norm is 
monotone for n — 1, but in general not for n > 2, as shown by Example 3 given in [9, 
Section 4]. A consequence of this example is that the proof of the nonexistence result 
for n — 1 cannot be extended to n > 2. 

In this paper we show that the nonexistence result itself cannot be extended. We 
show indeed that for n > 2 global-in-time classical solutions do exist, and also the 
persistence of supercritical regions fails. Our main result is the following (we state it 
for n — 2, but the generalization to any n > 2 is straightforward). 

Theorem 1.1 Let Q := {x <E M 2 : 1 < \x\ < 5} be an anulus. 

Then there exists u G C 2,1 (fi x [0, +oo)) satisfying equation (1.1) with T = +oo, the 
Neumann boundary condition (1.2), and 

{x G n : \Wu(x, 0)| > 1} = {x G D, : 2 < \x\ < 4} . 

Moreover there exists t > such that \Vu(x,t)\ < 1 for every (x,t) G Vl x (t , +oo). 

By the way, before Theorem 1.1 we didn't know any nontrivial example of (even 
local-in-time) transcritical classical solution in dimension n > 2. 

In order to prove Theorem 1.1 we can limit ourselves to radial solutions. Let r := \x\ 
be a radial variable, and let us consider radial solutions u(r, t). A simple computation 
(see [9]) shows that equation (1.1) becomes the following 

1 - ul 1 u r A , 

= WT^r"" + ;tt^' (l4) 

and the Neumann boundary condition (1.2) in the anulus becomes 

u r (l,t) = 0, (1.5) 

u r (5,f) = 0. (1.6) 
Therefore Theorem 1.1 is a consequence of the following result. 
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Theorem 1.2 There exists u G C 2,1 ([l,5] x [0, +00)) satisfying equation (1.4) in the 
strip [1, 5] x [0, +oo) ; the Neumann boundary conditions (1.5) and (1.6) for every t > 0, 
and the estimates 

0<w r (r,0)<l Vr G [1,2) U (4,5], (1.7) 

w r (r,0)>l VrG (2,4). (1.8) 

Moreover there exists t > such that \u r (x, t)\ < 1 for every (r, t) G [1, 5] x (t , +00). 

We recall that equation (1.4), without the second summand in the right-hand side, 
is just the Perona-Malik equation in dimension one. A bureaucratic count of derivatives 
says that the added term is a lower order term. Nevertheless its influence on the dy- 
namic is enormous. Due to that lower order term, the supercritical region of our radial 
solution disappears in a finite time, in contrast with the one dimensional case. After the 
extinction of its supercritical region, the solution becomes subcritical and has no more 
obstacles to global existence. For this reason what we need in the proof of Theorem 1.2 
is to keep the solution alive and regular up to this time. 

We are afraid that the possible extinction of supercritical regions makes the analytical 
study of the Perona-Malik equation even more difficult in dimension n > 2. We don't 
know whether this new phenomenon had been observed before in numerical experiments. 
We leave to numerical analysts any discussion about its consequences on the model and 
its practical applications. 

We conclude with some comments on our main result. 

Using an anulus instead of a ball should not be essential. This choice spares us and 
the reader from the technicalities due to the fact that equation (1.4) is singular for r = 0. 
The original equation however is not singular in the origin. Therefore this singularity 
should be easily compensated by the first Neumann boundary condition, which in a ball 
becomes u r (0, t) = 0. For this reason it should not be difficult to find solutions with the 
same properties, but defined in a ball. 

Concerning the nonlinearity, for simplicity we devoted this introduction to the model 
case of the Perona-Malik equation. On the contrary, in the following sections we work in 
a more general setting. What we actually do is to prove Theorem 1.1 and Theorem 1.2 
for the gradient flow of any integral functional with nonconvex integrand. We refer to 
the beginning of section 2 for the details. 

Finally, the reader could ask how special these solutions are. Theorem 1.1 indeed 
states the existence of just one such solution. What we actually construct is a large 
class of such solutions. In a word, they are not the fruit of some strange pathology, but 
a common feature in dimension n > 2. For example in the proof of Theorem 1.2, after 
choosing to small enough, we have some freedom in the choice of the initial condition 
u(r, 0) in the subcritical intervals [1, 2] and [4, 5], where we only impose some inequalities 
and some compatibility conditions at the endpoints (see section 3.1 for the details). Our 
construction then completes u(r, 0) in the supercritical interval [2, 4] in such a way that 
the solution starting with that datum is global and becomes subcritical for t > to- 
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This paper is organized as follows. In section 2 we reduce the construction of the 
required solution to the existence of solutions in four suitable subdomains, in each one 
of which the equation is either (degenerate) forward parabolic or (degenerate) backward 
parabolic. In section 3 we approximate the degenerate problems with strictly parabolic 
problems depending on a small parameter e > 0, and we state several e- independent 
estimates. In section 4 we prove these estimates. In section 5 we pass to the limit 
as £ — > + , showing that the limits are solutions of the degenerate problems. This 
completes the proof of Theorem 1.2, hence also of Theorem 1.1. 



2 The four subproblems 

Let us introduce some notations. Let T> C R 2 be a compact set. The parabolic interior 
of T> is the set Intp(P) of points (r, t) G T> for which there exists 5 > such that 
[r — 5, r + 5] x [t — 5, t] CD. The parabolic boundary of T> is the set dp(T>) : = T>\ Int p(T>). 

Throughout this paper we assume that (f G C°°(R) is an even function, hence in 
particular 

v/(0) = </'(0) = 0. (2.1) 

Moreover we assume that 

cp"(a)>0 Vae[0,l), (2.2) 

^(1) = 0, (2.3) 

ip"(a)<0 VaG(l,3], (2.4) 

(p'(3) > 0. (2.5) 

Of course the thresholds a = 1 and a = 3 can be replaced by any pair of positive 
numbers cr < o\. 

As a consequence of (2.2) and (2.3), or (2.3) and (2.4), we have that (f"'(l) < 0. As 
a consequence of (2.1) through (2.5) we have also that <p'(cr) > for every a G (0,3]. 
Figure 1 shows the typical behavior of <p'(o~). 




Figure 1: typical graph of ip' 
We consider the following equation 

u t (r,t) = y"{u r {r,t))u rr {r,t) + /Mr ' t)) , (2.6) 

r 
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which is the natural generalization of (1.4). 

It is easy to see that equation (2.6) reduces to (1.4) when <p(a) = 2 _1 log(l + a 2 ). 
We believe and we hope that this generality simplifies the presentation, and shows more 
clearly which properties of the nonlinear ity are essential in each step. 



2.1 The four regions 

In order to prove Theorem 1.2 we divide the strip [1, 5] x [0, +oo) into four regions. To 
begin with, we fix t > 0, and we consider the functions 



(3{t) := 3 - VI - t/to, 
defined for every t < to- Then we set 



7 (t) := 3 + v 7 ! - t/to 



Qi 

0,2 



0<t<t , 1 < r < /3(t)} , 
< t < t , P(t) <r< 7 (t)} 
< t < t , 7(t) < r < 5} , 



{(r,t) G 
{(r,t) G 
{(r,t) G 

[1,5] x [t ,+oo). 
We also set 

r a := {(/3(t),t) : t G [0,t ]} , T 3 := {{-y(t),t) : t G [0,t ]} . 

These sets are represented in Figure 2. We finally consider the functions 

^(1) 



b(t) :- 



mm 




xeR: ip"'{l)x 2 + /3'(t)x 







iftG [0,t ) 
if t = t , 



(2.7) 

(2.8) 
(2.9) 
(2.10) 
(2.11) 

(2.12) 
(2.13) 



c(t) . = , max ^xeR: V "\l)x 2 + 7 '(^ - ^ = j iftG[0,t ), (214) 

if t = t - 

It is not difficult to see that b(t) and c(t) are well defined for t G [0, to] provided that 
to is small enough. In Lemma 4.1 and Lemma 4.2 below we prove some properties of 
these functions, in particular their continuity. In Remark 2.5 we explain the heuristic 
idea behind these definitions. 

In order to prove Theorem 1.2 we need a solution of (2.6) in the strip [1, 5] x [0, +oo), 
satisfying the Neumann boundary conditions (1.5) and (1.6) for every t > 0, and es- 
timates (1.7) and (1.8). Our strategy is to construct this solution by glueing together 
solutions of the same equation in the four subdomains. These solutions are of course 
required to fulfil the Neumann boundary conditions (1.5) and (1.6). In order to be glued 
in a C 2 ' 1 way along the common boundaries, these solutions are also asked to satisfy 
the following conditions in F 1 

u(/3(t),t) = «(3, to) - 3 + (3{t) - <//(l) J ° -^y ds Vt G [0, t ], (2.15) 

u r ([3(t),t) = l VtG[0,t ], (2.16) 
u rr (P(t),t) = b(t) VtG[0,t ], (2.17) 

and the following conditions in T 3 

u(7(i),t)=«(3,fo)-3 + 7 (f)-^(l) f°^rds VtG[0,t ], (2.18) 

u r ( 7 (t),*) = l VtG[0,t ], (2.19) 

u rr (-f(t),t) = c(t) VtG[0,t ]. (2.20) 

Moreover we arrange things in such a way that equation (2.6) turns out to be (degen- 
erate) backward parabolic in Q2, and (degenerate) forward parabolic in the remaining 
three regions. Apparently all these conditions make the problem highly overdetermined. 
Nevertheless the degeneracy in Ti and T 3 allows us to fulfil all conditions. 



2.2 The four subproblems 

Let us state the existence results in the four regions. 

Theorem 2.1 (Region Qi) Let ip G C°°(R) be a function satisfying (2.1), (2.2), and 
(2.3). Let us assume that to > is small enough. Let j3(t), Qi, Ti, b(t) be defined by 
(2.7), (2.8), (2.12), and (2.13), respectively. 

Then there exists u G C 2,1 (<2i) satisfying equation (2.6) for every (r, t) G Q\, the 
Neumann boundary condition (1.5) for every t G [0,t ], and the boundary conditions 
(2.15) through (2.17). 

Moreover < u r (r,t) < 1 for every (r, t) G Qi \ IY 
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Theorem 2.2 (Region Q 2 ) Let ip G C°°(R) be a function satisfying (2.3), (2.4), and 
(2.5). Let us assume that t > is small enough. Let j3(t), ^{t), Q 2 , T l7 T 3 , b(t), c(t) 
be defined by (2.7), (2.9), (2.12), (2.13), and (2.14), respectively. 

Then there exists u G C 2,1 (Q,2) satisfying equation (2.6) for every (r,t) G Q.2 and 
the boundary conditions (2.15) through (2.20). 

Moreover u r (r, t) > 1 for every (r, t) £ Q 2 \ (I\ U r 3 ). 

Theorem 2.3 (Region <2 3 ) Let p G C°°(R) be a function satisfying (2.1), 
(2.3). Let us assume that t > is small enough. Let j(t), Q 3; T 3 , c(t) be 
(2.7), (2.10), (2.12), and (2.14), respectively. 

Then there exists u G C 2,1 (Q 3 ) satisfying equation (2.6) for every (r,t) 
Neumann boundary condition (1.6) for every t G [0, £o]> an d the boundary 
(2.18) through (2.20). 

Moreover < u r (r, t) < 1 for every (r, t) G <2 3 \ T 3 . 

Theorem 2.4 (Region Q 4 ) Let p G C°°(R) be a function satisfying (2.1), (2.2), and 
(2.3). Let t > be a real number, and let <2 4 be defined by (2.11). Let u G C 2 ([l,5]) 
be a function such that 

«or(l) = u 0r (5) = 0, (2.21) 
0<M 0r (r)<l Vr G [1,3) U (3,5]. (2.22) 

Then there exists a unique function u G C 2 ' l (Qi) satisfying equation (2.6) for every 
(r,t) G Qi, the Neumann boundary condition (1.5) and (1.6) for every t > t , and the 
initial condition u(r,t ) = uo(r) for every r G [1,5]. 

Moreover we have that < u r (r,t) < 1 for every (r,t) G Q4 \ {(3, to)}- 

The idea of constructing a solution by glueing solutions in suitable subdomains 
has already been used in [10]. The main difference is that in [10] all subdomains are 
rectangles (hence with fixed endpoints), and the prescribed values of u, u r , u rr at the 
boundary do not depend on time. 

Here the problem in <2 4 is a classical initial boundary value problem in a fixed 
interval. On the contrary, the problems in Q\, Q2, Q3 involve moving domains, no initial 
condition, but several time- dependent boundary conditions in the moving endpoints. 

As for initial conditions, in Q2 the backward character of the equation makes the 
solution completely determined by its values in Fi and T 3 . On the contrary, in Q 1 and 
<2 3 we have a lot of freedom in the choice of the initial condition, and for this reason 
there are plenty of different solutions. 

Concerning the multiple boundary conditions, let us consider for example the prob- 
lem in Q\. Once that an initial datum has been chosen, the Neumann boundary con- 
ditions (1.5) and (2.16) are enough to determine uniquely the solution. The surprising 
aspect is that this solution satisfies also (2.15) and (2.17), independently on the initial 
condition! In the following two remarks we show how the degeneracy of the equation in 
Ti makes this possible. 



(2.2), and 
defined by 

G Q 3; the 
conditions 
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Remark 2.5 Let u be a solution of equation (2.6) in Q 1 or Q 2 , with Neumann boundary 
condition (2.16). If u is of class C 2 ' 1 we can compute 



^K/3(t),t)] = / 9 , (tK(/3(t),t)+w t (/3(t),t) 

= ^K^) ) i)VKWt) ) t)W(t).t) + ^^ t) \ 

From condition (2.16) we have therefore that 

Integrating this equality in [t,t ] we obtain (2.15). This is actually a proof of (2.15). 
The proof of (2.18) is analogous. 

Remark 2.6 Let wbea solution of equation (2.6) in Qi or Q2, with Neumann boundary 
condition (2.16). Let us assume that u is smooth enough, and let us compute the time 
derivative of (2.16). We obtain that 

= j t [u r {(3{t),t)\ = f3'(t)u rr (P(t),t)+u rt (f3(t),t) 

2 ip"(u r )u rr (p'(u r ) 



= {3'(t)u rr + (p"(u r )u rrr + ip"'(u r )u rr + 



where in the last line all the derivatives of u are computed in ((3(t),t). Exploiting the 
Neumann boundary condition (2.16) we have therefore that 



P'(t)u rr (P(t),t) + ip"'(l)u 2 r Mt),t) - = 0. (2.23) 

This means that u rr (/3(t),t) is a solution of the equation defining b(t) in (2.13). If 
for t = we have that u rr ((3(0), 0) is the smallest solution of the equation, namely 6(0), 
then for all subsequent times u rr ((3(t),t) remains the smallest solution, namely b(t). The 
choice of the smallest solution is due to the fact that it tends to as t — > t , while the 
other solution diverges to +00. 

We point out that this simple argument is a heuristic justification of (2.17), but it 
is by no means a proof. In deriving (2.23) we used indeed that terms such as (p"(u r )u rrr 
vanish at the moving boundary. This requires some assumption on u rrr , which is beyond 
the C 2 ' 1 regularity. 
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2.3 Proof of Theorem 1.2 



Let us take to > small enough in order to apply Theorem 2.1, Theorem 2.2, and 
Theorem 2.3. Thus we obtain a solution of equation (2.6) in each one of the regions 
Qi, Qi, Qs- All these solutions are defined up to additive constants. We can therefore 
assume, without loss of generality, that they coincide in the common point (3, to). 

We claim that the three solutions glue together in order to give a solution of class 
C 2,1 in the rectangle [1,5] x [0,t ]. 

Let us examine indeed the solutions in Qi and Qi- If they coincide in (3, to), then 
they coincide in the whole r\ because they both satisfy (2.15). Also their first and 
second space derivatives coincide in Ti because they both satisfy (2.16) and (2.17). So 
these two solutions glue in a C 2,1 way. 

The same is true for the solutions in Q 2 and Q3. Note that in the common point 
(3, to) we have that u rr (3,to) = b(to) = c(t ) = 0. 

It remains to extend u to [1,5] x [t ,+oo). To this end, we apply Theorem 2.4 
using as "initial" condition the trace at t — to of the solution we have just glued in 
[1,5] x [0, to]- From Theorem 2.1 and Theorem 2.3 it is clear that this trace satisfies 
assumptions (2.21) and (2.22). This completes the proof of Theorem 1.2. □ 



3 Approximating problems 

In the previous section the proof of our main result has been reduced to the proof of 
Theorem 2.1 through Theorem 2.4. In this section we present our approach to Theo- 
rem 2.1 and Theorem 2.2. We skip Theorem 2.3 because it is symmetric to Theorem 2.1, 
and we skip Theorem 2.4 because it concerns a quite standard (non overdetermined) 
initial boundary value problem. 

The first thing to do is to choose to- Let us begin by considering the following 
constants depending only on <p 

7o :=3y9'(l) + 5, (3.1) 
71 := 5<//(l) + 100, (3.2) 
72 := max {\ip'(<j)\ + \<p"(a)\ + \<p"'(a)\ + \<p IV (<j)\ : a e [0,3]} . (3.3) 

Let us choose to € (0, 1) satisfying the following inequalities 

t0 -IRIjf' t0 -2^' t0 -96( 7l + l)V (3 ' 4) 

to ~ (20 7o 2 + 28 7o + 9)72' t0 ~ (12 7o +14) 72 - (3 ' 5) 

We stated these conditions on t as they are required throughout the proofs. It is 
not difficult to see that the third inequality in (3.4) implies the remaining four. 
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3.1 The forward problem in a moving domain 

Let us consider the problem in Q\. Let us choose uo G C°°([l,2]) satisfying the com- 
patibility conditions 

U0r(l) = 0, Uor(2) = l, u 0rr {2) = 6(0), (3.6) 

and the following inequalities 

0<M r(r)<l We (1,2), (3.7) 

Kr(r)|<10 We [1,2], (3.8) 

|w 0rrr (r)| < 10 We [1,2]. (3.9) 

We remind that 6(0) depends on to, and for this reason the choice of uo depends on 
the choice of to. 

Now we approximate the problem in Q 1 with a strictly parabolic problem depending 
on a small parameter e. 

Theorem 3.1 Let if G C°°(R) be a function satisfying (2.1), (2.2), (2.3). Let us 
assume that t > satisfies (3.4) and (3.5). Let j3(t), Q 1; b(t), be defined by (2.7), 

(2.8) , and (2.13), respectively. Let u G C°°([l,2]) be a function satisfying (3.6) through 

(3.9) , and lets G (0,1). 

Then there exists a unique function u e G C 2,1 (Qi) satisfying equation (2.6) for every 
(r,t) G Qi, the Neumann boundary condition (1.5) for every t e [0,to]; the further 
Neumann boundary condition 

u £ r (P(t),t) = l-e Vfe[0,* ], (3.10) 

and the initial condition 

u £ (r,0) = (1 -e)uo(r) Vr G [1,2]. 

Moreover u £ satisfies the following estimates independent on e. 

1. Maximum principle for space derivatives. We have that 

<u £ r (r,t) <l-e V(r,t)GQi. (3.11) 

2. Uniform strict parabolicity in the interior. For every 5 G (0, 1) there exists a 
constant Mi = Mi (5) such that 

u £ (r, t) < Mi < 1 V£ G [0, t ], Vr G [1, f3{t) - 5]. (3.12) 

As a consequence, there exists a constant M 2 = M 2 (S) > such that 

p"(u £ (r,t)) >M 2 \/t G [0,t ], Vr G [l,(3(t) - 5}. (3.13) 
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3. Estimate on second derivatives at the fixed boundary. We have that 

< u e „(l,t) < 100 Vfe[0,f o ]. (3.14) 

4- Estimate on second derivatives at the moving boundary. We have that 

\u £ rr (P(t),t)-b(t)\<e Vfe[0,* o ]. (3.15) 

5. Global estimate on second derivatives. There exists a constant M 3 such that 

\u £ rr (r,t)-b(t)\<e + M 3 \r-f3(t)\ V(r,t)eQi- (3.16) 

As a consequence, there exist constants M 4 and M 5 such that 

\u £ rr (r,t)\<M 4 V(r,*)eQi, (3.17) 

\u e t (r,t)\ < M 5 V(r,*)eQi- (3.18) 

6. Integral estimates. There exists a constant M 6 such that 

<p"(u £ r (r,t))[u £ rt (r,t)} 2 dr dt < M 6 . (3.19) 



/ 

Jo. 



In addition, for every 5 G (0, 1) there exists a constant M 7 = M 7 (5) such that 
rf3(t)S 

J [u £ rt {r,t)f dr < M 7 VtG[0,t ], (3.20) 

J [u £ rrr (r,t)} 2 dr < M 7 VtG[0,t ], (3.21) 

t r /3(t)-5 

[u £ rrt (r,t)} 2 drdt< M 7 . (3.22) 



3.2 The backward problem in a moving domain 

The backward problem in Q 2 can be transformed in a forward parabolic problem by 
reversing the time. To this end we consider the domain 

T := {(r,t) G R 2 : < t < t , f3(t - t) < r < t(*o -*)} • (3.23) 

The domain T is just Q 2 upside-down, while its parabolic boundary dpT is just 
T 1 U T 3 upside-down (see Figure 3). 

With this notation proving Theorem 2.2 is equivalent to proving the following result 
(note in particular the minus signs in the right-hand side of (3.24)). 
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Figure 3: domains Q 2 , T, and % 



Theorem 3.2 (Region T) Let <p £ C°°(M) be a function satisfying (2.3), (2.4), and 
(2.5). Let us assume that t > satisfies (3.4) and (3.5). Let Pit), jit), T, b(t), c(t) 
be defined by (2.7), (3.23), (2.13), and (2.14), respectively. 
Then there exists u £ C 2,1 (T) satisfying equation 

u t (r,t) = -p''MM))urrM) - ^ (Mr(r ' t)} (3.24) 

r 

for every (r,t) £ Q2, and the Neumann boundary conditions 

u r ((3(t -t),t) = u r ( 7 (t -t),t) = l Vt £ [0, t \. 

Moreover we have that 

u r (r,t)>l V(r,t) £ Intp(T), 

U rr (f3(to-t),t)=b(t -t) Vt £ [0,* ], 

Mrr ( 7 (t -t),t) = c(t -t) Vt£ [0,t ]. 

We didn't state explicitly what (2.15) and (2.18) become in this new setting because, 
as we have seen in Remark 2.5, they easily follow from the regularity of the solution and 
the Neumann boundary conditions (2.16) and (2.19). 

In order to approach Theorem 3.2, for every e £ (0, to) we se t (see Figure 3) 

% := {(r,t) £ T : t > e} . (3.25) 

Now we approximate the degenerate parabolic problem in T with a strictly parabolic 
problem in T £ . Note that in this case we prescribe an initial condition for t = e. This 
initial condition is compatible with the Neumann boundary conditions. 

Theorem 3.3 Let ip, t , (3{t), j(t), T, b(t), c(t) be as in Theorem 3.2. Lets £ (0,t ), 
and let % be defined by (3.25). 

Then there exists a unique function u £ £ C 2 ' l {T £ ) satisfying equation (3.24) for every 
(r,t) £ T e , the Neumann boundary conditions 

u £ r ((3(t -t),t)=u £ Mt -t),t) = l + e Vt £ Mo], 
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and the initial condition 

u e (r, e) = (1 + e)r Vr G [/3(t - e), 7(*o - e)]- 

Moreover u £ satisfies the following estimates independent on e. 

1. Maximum principle for space derivatives. We have that 

1 + e < u £ (r, t) < 3 V(r, i) G %. (3.26) 

2. Uniform strict parabolicity in the interior. For every compact set K C T^fllntp(T) 
there exists a constant M x ( depending on K but not on e) such that 

1< Mi < <(r, t) V(r, t) G X. (3.27) 

As a consequence, there exists a constant M 2 > such that 

- <p"(u e r (r, t)) > M 2 V(r, t) G X. (3.28) 

3. Estimate on second derivatives at the boundary. We have that 

\u £ rr ((3(t -t),t)-b(t -t)\<^ Vfe[e,to], (3.29) 
KMto-t),t)-c(t -t)\<y/i Wte[e,t ]. (3.30) 

4- Global estimate on second derivatives. There exists a constant M 3 such that 

\u £ rr (r,t)-b(t -t)\ <y/£ + M 3 \r-P(t -t)\ V(r,t)e%, (3.31) 

\u e „(r,t)-c(t -t)\< y /e + M 3 \r-'y(t -t)\ V(r,t)e%. (3.32) 

As a consequence, there exist constants M4 and M 5 such that 

\u £ rr {r,t)\ < M 4 V(r,t)eT e , 
\u £ t (r,t)\ < M 5 V(r,t)e%. 

5. Integral estimates. There exists a constant M 6 such that 

- [ ip"(u £ r (r,t))[u £ t (r,t)} 2 drdt< M 6 . (3.33) 

In addition, for every closed rectangle [ri,r 2 ] x t 2 ] ^ ^ H Intp(T) i/iere exists 
a constant M 7 ( depending on the rectangle, but not on e) such that 

rr-2 

/ [u £ rt (r,t)] 2 dr <M 7 Vte[t u t 2 ], 
J n 
/•f 2 

/ [u £ rrr (r,t)] 2 dr < M 7 Vte[ti,t 2 ], 

J ri 

I' [ 2 [u £ rrt (r,t)} 2 drdt<M 7 . 
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4 Proof of estimates 

In the following two statements we collect the properties of b(t) and c(t) which are 
needed in the main proofs. We only prove Lemma 4.1 because the proof of Lemma 4.2 
is quite similar. 

Lemma 4.1 Let t > 0, and let (5{t) be defined by (2.7). Let us assume that 

f > 4^(1)1^(1)1. (4.1) 
to 

Then the function b(t) introduced in (2.13) is well defined for every t G [0, t ] and 
fulfils the following properties. 

(1) We have that b G C°([0, t }) n C°°([0, t )). 

(2) We have that 

< -b'{t) < 5t cp'(l)f3\t) V* G [0,t ). (4.2) 

(3) The function b{t) is decreasing and 

< b(t) < b(0) < 2<p'(l)t VtG[0,t ], (4.3) 

< 6(* - t) < 2<p'(l)y/hVi Vte[0,* o ]. (4.4) 
Proof. From (2.7) we easily deduce that 

P(t) > 2, (3'(t) > (3"(t) = 2t [P\t)f (4.5) 

for every t G [0, to]. Moreover from assumption (4.1) we have that 



[/n*)] 2 >^>4yy(l)|yAl)l VtG[0,t ), (4.6) 



Let us consider the polynomial 



p(x) :=if" / (l)x 2 + (3'(t)x 



where t G [0,t ] is thought as a parameter. Due to our assumptions on ip we have that 
(p"'(l) < 0. Exploiting (4.6) one can therefore check that 
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This implies that p(x) has at least one real root. Moreover p(x) is either a polynomial 
of degree one, or a polynomial of degree two representing a concave function. In both 
cases the unique or the smallest root b(t) is the one lying in the interval whose endpoints 
appear in (4.7), hence 

^ <b(t)<^- Vte[0,t ). (4.8) 



(3\t)(3'{t) " w " (3>(t) 

The regularity of b(t) in [0, t ) follows from the regularity of the coefficients of the 
equation. The continuity up to t — t follows from (4.8) and the fact that (3'(t) — > +oo 
as t — > (we remind that we set b(t ) = 0). This proves statement (1). 

Let us consider now the derivative b'(t). From the implicit function theorem we have 
that 

Let us estimate the numerator of (4.9). Exploiting (4.5) and the upper bound for 
b(t) provided by (4.8), we have that 

o < w(o+2^ja < 2 tom)] ^ + qi m 



4/3'(t) 

Let us estimate the denominator of (4.9). Exploiting once again the upper bound 
for b(t) provided by (4.8), and inequality (4.6), we obtain that 

2^(1)6(0 + m > m - 2i^(i)i^|| = m (i - 2| ^j}p (1) ) > \m. 

At this point (4.2) easily follows from the estimates on the numerator and the de- 
nominator of (4.9). This completes the proof of statement (2). 

Let us prove now statement (3). The function b(t) is decreasing because its derivative 
is negative. It follows that = b(t ) < b(t) < 6(0) for every t G [0, to]- Finally, the 
estimates on 6(0) and b(t — t ) stated in (4.3) and (4.4) follow from (4.5) and the upper 
bound in (4.8). □ 

Lemma 4.2 Let to be as in Lemma 4-1, and let ^(t) be defined by (2.7). 

Then the function c(t) introduced in (2.14) is well defined for every t e [0, to] an d 
fulfils the following properties. 

(1) We have that c e C°([0, t ]) n C°°([0, to)). 
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(2) We have that (we remind that 7'(t) is negative) 

< c'(t) < -5t ¥>'(l)Y(t) 



Vt e [0,t ). 



(3) The function c(t) is increasing and 

- 2ip'(l)t < c(0) < c(t) < 
> c(t - 1) > -2cp'(l)VhVt 



Vt G [0, t ], 
Vt G [0,t ] 



(4.10) 
(4.11) 



Remark 4.3 The second inequality in (3.4) implies (4.1). We can therefore apply the 
conclusions of Lemma 4.1 and Lemma 4.2 in the proofs of our main results. Exploiting 
the first and second inequality in (3.4), from statement (3) in Lemma 4.1 and Lemma 4.2 
we obtain the following stronger inequalities 



We finally state a classical comparison result for fully nonlinear parabolic equations. 
This is the key tool in our analysis. We omit the standard proof. 

Lemma 4.4 Let V C M 2 be a compact set, and let ip : V x IR 3 — > R be a continuous 
function. Let us assume that 

• if) is nondecreasing in the last variable (degenerate ellipticity), namely 

i/j(r, t, p, q, Si) < ij}{r, t, p, q, s 2 ) Vsi < s 2 , V(r, t,p,q) E V x R 2 , 

• tp is locally Lipschitz continuous in the third variable, namely for every R > 
there exists a constant L such that 

\tp(r,t,pi,q,s) -*/>(r,t,p2,q,s)\ < L\p x -p 2 \ V(r, t,p u p 2 , q, s) eVx [0,R] 4 . 
Let u and v be two functions in C°(V) n C 2,1 (lntp(V)) such that 



Then u(r,t) < v(r,t) for every (r,t) G V. □ 

Remark 4.5 If we know a priori that there exist constants a and b such that a < 
u(r,t) < b for every (r, t) G T>, then we can weaken the assumptions of Lemma 4.4 by 
asking that v fulfils its differential inequality only for those (r, t) G Intp('D) such that 
a < v(r, t) < b. We can obtain a similar statement by swapping the role of u and v. 



< b{t) < 6(0) < 1 
< b(t -t)<y/t 
-Vt < c(t - t) < 



Vt G [0, t ], 
Vt G [0, t ], 
Vt G [0, t ] 



(4.12) 
(4.13) 
(4.14) 



u t < ip(r, t, u, u r , u TT ) 

v t > ^(r,t,v,v r ,v rr ) 
u(r, t) < v(r, t) 



V(r,t) Glntp(P), 

V(r,t) G Intp(P), 
V(r,t) ed P (V). 
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4.1 Proof of Theorem 3.1 

Let us briefly sketch the outline of the proof. First of all we show that the initial 
boundary value problem has a unique solution u £ G C 2,1 (Qi) fl C°°(Intp(Qi)). Then 
for the sake of simplicity we set v =: u £ r , and w := u e rr . These functions belong to 
C°(<2i) fl C°°(Intp(Qi)). It is easy to see that v is the solution of the following equation 

v t = tp"(v)v rr + <p"'{v)v 2 r + ^-v r - ^ V(r,t) G Int P (Q 1 ), (4-15) 
with Dirichlet boundary conditions 

v(l,t) = Vte[0,f o ], (4.16) 

v(P(t),t) = l-e Vfe[0,t ], (4.17) 

and initial datum 

u(r,0) = (1 -e)uor(r) Vr G [1,2]. (4.18) 
In the same way w is a solution in Intp(Qi) of equation 

w t = (p"(v)w rr + 3(f'"(v )w r w + (p IV (v)w 3 

ip"'(v) 2 <p"(v) <p"(v) <p'(v) /A s 

/y rjrt /y A /yO 

where the terms in v are to be interpreted as coefficients depending on r and t. Moreover 
w satisfies the Dirichlet boundary conditions 

w(l,t) =< r (l,t) Vte[0,t ], (4.20) 

w((3(t),t) = u £ rr ((3(t),t) Vt G [0,t ], (4.21) 
and the initial condition 

w (r, 0) = (1 - e)u Qrr {r) Vr G [1, 2]. (4.22) 

The initial data (4.18) and (4.22) can be easily estimated using (3.6) through (3.9). 
Indeed from Taylor's expansion we have that 

u 0r {r) = u 0r (2) + u 0rr {2)(r - 2) + ^u 0rr r(0( r ~ 2 Y 
for a suitable £ G (1,2), hence 
1 + 6(0)(r - 2) - 5(r - 2) 2 < u 0r (r) < 1 + 6(0)(r - 2) + 5(r - 2) 2 Vr G [1, 2]. (4.23) 
Analogously we have that uo rr (r) = uq„(2) + Morrr(0( r — 2), hence 

6(0) - 10(2- r) < u 0rr .(r) < 6(0) + 10(2- r) Vr G [1,2]. (4.24) 

Then we show that t> and w satisfy five sets of inequalities in Q\. 



17 



• The first pair of inequalities is 

< v(r,t) < 1 -e, (4.25) 

which is equivalent to (3.11). 

• The second inequality is 

v(r,t) < l-7]^(r-3) 2 + ^-l^ (4.26) 

for a suitable constants rj > 0. The term after rj is positive in Q^IY This implies 
(3.12), hence also (3.13). 

• The third pair of inequalities is 

2D 

^<^ (l-800 7 ^ (l - e (4 - 27) 

where 72 is the constant defined in (3.3). The left-hand side and the right-hand 
side coincide for r = 1, hence v r (l,t) is bounded by their derivative computed in 
r = 1. Thus we obtain that < u r (l,t) < 20(1 - 8007 2 t)~ 1/2 . The upper bound 
is less than 100 due to the second inequality in (3.4). This proves (3.14). 

• The fourth pair of inequalities, and probably the most delicate one, is 

v(r, t)>l-e- (b(t) + e)(J3(t) - r) - 7 o(/3(t) - rf, (4.28) 

v(r, t)<l-e- (b(t) - e){/3(t) - r) + lo (/3(t) - r) 2 , (4.29) 

where 70 is the constant defined in (3.1). In this case v(r,t) is bounded from 
below and from above by two functions which coincide for r = (3(t). It follows 
that v r (/3(t),t) is bounded by the space derivatives of these two functions computed 
in r = /3(t). This implies (3.15). 

• Thanks to (3.14) and (3.15) we have an estimate on the values of w at the bound- 
ary. This is the starting point to prove the fifth set of inequalities 

6(f) - e - 7 i(/3(t) - r) < w(r, t) < b(t) +e + 7 i(/3(i) - r) (4.30) 

where 71 is the constant defined in (3.2). This implies (3.16), hence also (3.17) 
and (3.18). 

All these estimates are proved using subsolutions and supersolutions. Finally the 
proof of (3.19) through (3.22) relies on usual energy estimates. 
We are now ready to proceed with the details. 
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Existence and maximum principle for space derivatives Let us take a function 
(f £ G C°°(R) which coincides with (p in the interval [0, 1 — e\, and such that (p"(cr) > 
v e > for every trel Equation (2.6), with ip E instead of </?, is strictly parabolic. The 
initial boundary value problem in Qi can be reduced to the fixed domain [0, 1] x [0,to] 
by the variable change (s,t) — > (1 + (Pit) — l)s,t), which of course doesn't change the 
strict parabolicity. By well knows results (see for example [16]) the problem admits a 
unique solution u e G C 2,1 (<2i) fl C 00 (Intp(Qi)). By the way, the estimates we are going 
to prove could be used to give a self contained proof of the existence result for e fixed, 
showing in particular that the singularity of Q\ in (3, to) doesn't affect the existence or 
the regularity of the solution. 

If we show that this solution satisfies (3.11), then this same solution satisfies equation 
(2.6) with the original ip. Our assumptions on <p imply that <p'(0) = and ip'(l — e) > 0. 
It follows that z(r,t) := is a subsolution of (4.15) through (4.18), while z(r, t) := 1 — £ 
is a supersolution of the same problem. The usual comparison principle implies (4.25). 

Uniform strict parabolicity in the interior Let us choose rj > small enough so 
that 

Note that the infimum is positive because our assumptions on u imply that the 
fraction is positive for every r£ [1,2) and tends to 6(0) /2 > as r — > 2~. 

Let z(r,t) denote the right-hand side of (4.26). We claim that, for this choice of r], 
z is a supersolution of (4.15) through (4.18), which implies (4.26). 

Boundary r = 1 Since \r — 3| < 2, t < t , and rj < 1/8, we have that 

z(r, t) = 1 - V (V - 3) 2 + ^ - 1^ > 1 - 4 V > l - V(r, t) e Qi- (4.32) 
In particular we have that z(l,t) > = v(l,t). 

Boundary r = /3(t) In this case we have that z(/3(t),t) — 1 > 1 — e — v(/3(t),t). 

Boundary t = Thanks to the second inequality in (4.31), and the fact that 
Uor(r) > for every r G [1, 2], we have that 

z( r , 0) = 1 - 7] ((r - 3) 2 - l) = 1 - r]{2 - r)(4 - r) > 

> u 0r {r) > (1 - e)u 0r {r) = v(r, 0). 
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Differential inequality We have that 

zt = -j-, z r {r, t) = -2i](r - 3), z rr (r, t) = -2rj. 
to 

Therefore we have to prove that 

-- > -2 W "{z) + <p"'(z)4ri 2 (r - 3) 2 - ^-2 V (r - 3) - 
to r r 2 

Due to (4.32) and the properties of tpf , we have that (p'(z)/r 2 > (l/9)<//(l/2). All 
the other terms are uniformly small when rj is small. This shows that the differential 
inequality is satisfied when rj is small enough, for example as soon as i] fulfils the last 
inequality in (4.31). 

Estimate on second derivatives at the fixed boundary Let us prove (4.27). We 
already know from (4.25) that v(r,t) > in Q-y. In order to prove the other inequality, 
we set for simplicity k(t) := 20(1 — 8007 2 t)~ 1 / 2 . It is easy to check that k(t) is the 
solution of the Cauchy problem 

k'(t)= l2 k\t), fc(0) = 20. 

We claim that z(r,t) := k(t) (1 — e 1_r ) is a supersolution of (4.15) through (4.18), 
which implies (4.27). 

Boundary r = 1 In this case we have that z(l,t) = = v(l, t). 

Boundary r = f3(t) Since f3(t) > 2, in this case we have that 

z((3(t),t) = k(t) (1 - e 1 -^) > 20(1 - e' 1 ) > 1 > v((3(t),t). 

Boundary t = From the assumptions on uq we have that 

v(r, 0) = (1 - e)u 0r (r) < u 0r (r) = u 0r (l) + (r - l)«orr(0 < 10(r - 1). (4.33) 
On the other hand, the function 1 — e 1_r is concave, hence 

z(r, 0) = 20 (1 - e 1 ^) > 20(1 - e" 1 )^ - 1) > 10(r - 1). (4.34) 
From (4.33) and (4.34) it follows that z(r, 0) > v(r,0) for every r e [1,2]. 
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Differential inequality We have that 

Zt (r, t) = k'(t) (1 - e 1 -") , z r (r, t) = k^e 1 -' , z rr {r, t) = -A;(t)e 1 " r . 
We have therefore to prove that 

fc'(t) (1 - e l ~ r ) > V "'(z)k 2 (t)e 2 - 2r - k(t) (l - ^ (p"(z)e 1 - r - (4.35) 

Thanks to Remark 4.5, this inequality has to be satisfied only for those values (z, t) G 
Qi for which < z(r, t) < 1. In this case we know that the last two terms are negative. 
Since <£>"'( 0) = we have also that 

W"{z)\ = W"{z) - v/'(0)| = z < l2 k(t) (1 - e 1 -) . 

It follows that 

right-hand side of (4.35) < \^'"{z)\ k 2 {t) < l2 k 3 {t) (l - e^') = k\t) (l - e 1 ^) , 
which completes the proof of the differential inequality. 

Estimate on second derivatives at the moving boundary Subsolution Let 

us prove (4.28). To this end, we denote its right-hand side by z(r,t), and we show that 
z is a subsolution of (4.15) through (4.18). 

Boundary r = 1 Since (3(t) > 2 and 7 > 1, in this case we have that 

z(l,t) = l-e- (b(t) + e)((3(t) - 1) - 7o (/3(f) - l) 2 < 1 - 7 o < = v(l, t). 

Boundary r — (5{t) In this case we have that z((3(t),t) = 1 — e — v(/3(t), t). 

Boundary t = Since 6(0) > and 7 > 5, from (4.23) have that 

v(r, 0) = (1 - e)u 0r (r) > (1 - e) (l - 6(0)(2 - r) - 5(2 - r) 2 ) > 
> 1 - e - (6(0) + e)(2 - r) - 7 (2 - r) 2 = z(r, 0). 

Differential inequality Let us set for simplicity x := (3{t) — r. Then we have that 
z r (r, t) = b(t) + e + 2^ x, z rr (r, t) = -2~f , 

z t (r,t) = -b(t)P'(t) - (3\t)e - (b'(t) + 2 l0 (3\t))x. 
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The differential inequality is therefore the following 



-b(i)0(i) - 0(i)e - (b'(t) + 2 l0 p'(t))x < -2 w "(z) + f'"(z)(b(t) + e + 2 l0 xf 



Thanks to Remark 4.5 we can limit ourselves to prove it for all values (z, t) e Q\ 
such that < z(r, t) < 1. From the definition of b(t) we have that 



b{t)(3\t)=^"\l)b\t) 



(4.36) 



After changing the signs we have therefore to prove that 

P'(t)e + (b'(t) + 2 l0 p'(t))x > 2 lQ <p'\z) 

-y\ z )(b{t)+e + 2 lG xf-^"{l)b\t)\ 
ip"{z) 



+ 



■(b(t)+e + 2 lo x) 



<p'(z) <//(!) 



- 2 p 2 (t) 
h + I 2 + h + h- 



(4.37) 



Let us estimate the left-hand side. Exploiting (4.2), and the fact that t$ < I and 
27o > 2 + 5</?'(l), we have that 



f3'(t)e + (b'(t) + 2 lo (3'{t))x > f3'(t)e + (-5<p'(l)t + 2 7o ) (3'{t)x 

> p'(t)e + 2p'(t)x 



^ £ X 

2to to 



(4.38) 



Let us estimate the four terms in the right-hand side of (4.37). From now on we 
exploit several times the following inequalities 

< e < 1, < x < 2, < z < 1, < b(t) < 1, 1 < r < (3{t). 

First of all we remark that 

\z-l\ <e + {b{t) + e + 7o x)x < e + (2 7o + 2)x. 

As a consequence we have that 

\<p"(z)\ = \ip"(z)-ip"(l)\ = \z-l\- |v? w (OI < 126 + 72 (2 7 o + 2)x, (4.39) 

hence 



h < 2 7o |^"(z)| < 2 727 oe + 72(47^ + 4 7o )x. 



(4.40) 
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Analogously we have that 

\<ff"(z) - y/'(l)| = \z - 1| • |/ V (0I < 126 + 72 (2 7o + 2)x, 



and 



hence 



| (6(f) + e + 2 lo x) 2 - b 2 (t) | = (26(f) + e)e + (4 7o 2 a; + 4&(t) 7o + 4e 7o )x 

< 3e + (8 7o 2 + 8 7o )x, 



h < \<p"'(z)\-\(b(t)+e + 2 lo x) 2 -b 2 (t)\+b 2 (t)\<p"'(z)-<p"'(l)\ 
< 4 72 e + 72 (8 7 o + 10 7o + 2)x. 



Exploiting (4.39) once again, we obtain that 



h < (b(t)+e + 2 lo x) 



< (47o + 2)|^)| 

< 72 (4 7o + 2)e + 72 (8 7 ^ + 12 7o + 4)x. 

Finally we have that 

\<p'(z) - </?'(!) | = \z-l\- \tp"(0\ < lie + 72 (2 7o + 2)x, 



and 



hence 



1 1 



P 2 (t) 



\(3(t)+r\-\(3(t)-r\ . \m+r\ , 

„ „ , , Ob Ob ^ 



h < W(z)\ 



r 2 f3 2 (t) 



1 1 



- R2 



r 2 f3 2 (t) 
< i2£ + 72 (2 7o + 3)x. 



P 2 (t) 



f3 2 (t) 



\</(z)-</(l)\ 



Summing (4.40) through (4.43), and exploiting (3.5), we obtain that 



e x 



2to to 



This estimate and (4.38) imply the differential inequality. 



(4.41) 



(4.42) 



(4.43) 



h + h + h + h< 72 (6 7o + 7)e + 72 (20 7o 2 + 28 7o + 9)x < — + -. (4.44) 



Estimate on second derivatives at the moving boundary - Supersolution Let 

us prove (4.29). To this end we show that its right-hand side, which we denote by z(r, t), 
is a supersolution of (4.15) through (4.18). 
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Boundary r = 1 Since 70 > 1 > b(t) and f3(t) > 2, in this case we have that 
z(l, t) = l-e+ (P(t) - 1) [ 70 (/3(0 - 1) - b(t) + e] > 1 - e = v(l, t). 

Boundary r — j3{t) As in the case of the subsolution we have that 

z((3(t),t) = l-e = v((3(t),t). 

Boundary t = Since 70 > 5 and 6(0) < 1, from (4.23) we have that 
v(r, 0) = (1 - e)u 0r (r) < (1 - e) [l + 6(0)(r - 2) + 5(r - 2) 2 ] < 
< 1 - e + (6(0) - e)(r - 2) + 7o(r - 2) 2 = z(r, 0). 

Differential inequality We limit ourselves to sketch the argument, which is analo- 
gous to the case of the subsolution. After computing the derivatives and using equation 
(4.36), we reduce ourselves to prove that (as before x := f3(t) — r) 

P , (t)e + (2 l0 p'(t)-b'(t))x > 2 W "(z) 

+ y\ z )(b{t)-e-2 lQ xf-^"{l)b\t)] 



+ 



-(&(t)-e-2 7o a;) 

^(1) <p'(z) 
(3 2 (t) r 2 



=: h + I 2 + h + h. 
Let us estimate the left-hand side. Since b'(t) < and 70 > 1, we have that 

P'(t)e + (270/7(0 - b'(t))x > P'(t)e + 2p(t)x > ^ + f . 

AIq to 

The four terms in the right-hand side can be estimated as in the case of the subso- 
lution. In this way we obtain (4.44) once again. We conclude by using the smallness of 
to exactly as in the case of the subsolution. 



Global estimate on second derivatives — Supersolution Let us prove the upper 
bound in (4.30). To this end we denote the right-hand side by z(r, t), and we show that 
it is a supersolution of (4.19) through (4.22). 

Boundary r = 1 Since b(t) > 0, f3(t) > 2, and 71 > 100, from (3.14) we have that 
z(l,t) = b(t)+e + 7l (/3(t) - 1) > 71 > 100 > u £ rr (l,t) = w(l,t). 
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Boundary r = j3(t) From (3.15) we have that 

z(P(t),t) = b(t)+e> u e „(P(t),t) = w(P(t),t). 

Boundary t = Since 7! > 10, from (4.24) we have that 

z(r, 0) > 6(0) + 7i (2 - r) > 6(0) + 10(2 - r) > (1 - e) [6(0) + 10(2 - r)] > 

> (1 - e)u 0rr (r) = w(r, 0). 

Differential inequality We have that 

z r (r, t) = -7i, z rr {r, t) = 0, z t {r, t) = b'(t) + ^f3'(t). 
The inequality to be satisfied in Qi is therefore the following 

b'it) + 7l /7(f) > -^"'(v)z + ^(,), 3 + ^M, 2 - 71^ - + 2^M, 

where of course the terms in v are thought as coefficients. 

Let us examine the left-hand side. Using (4.2) and the fact that t < 1 and 71 > 
5<^'(1) + 1, we have that 

b\t) + 7l /3'(0 > (-5^'(l)t + 7i) P(t) > (3'(t) > ± (4.45) 

Since \z(r, t)\ < 2 + 271, a rough estimate of the right-hand side of the differential 
inequality gives that 

right-hand side < 72(471 + 6) (2 + 27i) 3 < 4872(71 + l) 4 . 

Due to the last condition in (3.4), this estimate and (4.45) imply the differential 
inequality. 

Global estimate on second derivatives — Subsolution Let us prove the lower 
bound in (4.30). To this end we denote the left-hand side by z(r, t), and we show that 
it is a subsolution of (4.19) through (4.22). 

Boundary r = 1 Since (3(t) — 1 > 1 and b(i) < 1 < 71, from (3.14) we have that 
z(l, t) = bit) - e - 7 i(/3(t) - 1) < 1 - 71 < < u e „{\, t) = w{\, t). 

Boundary r = (3(t) From (3.15) we have that 

z((3(t),t) = b(t) - e < u e „(0(t),t) = w(P(t),t). 
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Boundary t = Since 6(0) < 1 and 71 > 10, from (4.24) we have that 
z( r , 0) = 6(0) - e - 71(2 - r) < (1 - e)(6(0) - 10(2 - r)) < (1 - e)u 0rr {r) = w(r, 0). 

Differential inequality We have to prove that 

b'tt) - - n m < w»* + jy {v y + + 7l ^l - + 2^i. 

Since 71 > 1 and < 0, in the left-hand side we have that 

- 7i/7(f) <-/?(*)< 
while the usual rough estimate on the right-hand side gives that 

right-hand side > —4872(71 + l) 4 . 
The conclusion follows as in the preceding case. 

Integral estimates In the following estimates we introduce constants ci, C2, . . . , all 
independent on e. 

Computing the time derivative and integrating by parts, we have that 



dt 



I [utfdrj = 2^ ul--(r V '(ul)) rt dr + (3'(t)[umt),t)] 2 

= -2 J </«) K,] 2 dr + 2j ^f±u\u\ t dr 



+2<p"(l-e)ut(f3(t),tKMt),t)+P'(t) HW),t)] 2 
=: h + I 2 + h + h. (4.46) 

Note that in the integration by parts we neglected the boundary term in r = 1 
because u £ rt (l,t) = due to the Neumann boundary condition (1.5). 

Let us estimate some of the terms in (4.46). From (3.18) we have that 

< </«) ([<f + K,] 2 ) < y/«) [< t f + Cl , 

hence 

rP{t) 

h< J v"K)Kt] dr + c 2 . (4.47) 
From (3.18) we have also that 

h < c 3 p'(t). (4.48) 
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Taking the time derivative of the Neumann boundary condition (3.10) we obtain 
that 

From (3.17) we have therefore that 

\uUm,t)\=P'(tMMt),t)\<c 4 (3 f (t), 

hence by (3.18) 

h < c 5 P'(t). (4.49) 
Plugging (4.47), (4.48), and (4.49) in (4.46) we obtain that 

^ ( jf K] dr ]^-J i ^K) Kt\ dr + c 2 + c 6 p'(t), 

and therefore 

rto rP(t) z-2 

/ / <p"(u £ r )[u £ rt } 2 drdt< K(r,0)] 2 rfr + c 2 t + c 6 (/3(to)-/5(0))<c 7 . 
Jo Ji Ji 

This proves (3.19). 

In order to prove (3.20) through (3.22), we choose a cut-off function p e C°°(R 2 ) 
such that < p(r,t) < 1 for every (x,t) G IR 2 , and 

p(r,t) = 1 in{(r,t)eQi:l<r</3(t)-<J}, 
p(r,t) = in{(r,OeQi:/3(0-<V2<r</3(0}- 

Then we set 

rW) f P(t) 
E(t):=J p 2 K t ] 2 dr, F ^ : =J 1 P 2 ^'\<)[<rtfdr, 

G(t):= ^ V "(u £ r )[u £ rt ] 2 dr. 
Taking the time derivative of E(t) and integrating by parts we have that 

,/3(t) r l3(t) 



pp t [u £ t ] 2 dr + 2 J p 2 u £ rt u e rtt dr 

,/?(t) ,/?(t) 
= 2 J PPt[u £ rt } 2 dr - 2 y (p 2 u £ rt ) r u e tt dr, 
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where we neglected the boundary terms because u rt = when r = 1, and p = when 
r = P(t). Computing the derivatives in the right-hand side we end up with 

,/3(t) 

dr 



E'(t) = 2 J p Pt [u £ J 

r p(t) f p(t) 
-4 J p PrV '"(K) dr-Aj pp rV "{ul)ul t ul rt dr 

-4 J pPr 1S-lly u ^ dr-2j pV'KK t ulX r , dr 
-2 J pVK) K r ,f dr-2j^ p 2 ^-uUK rt dr 



=: h(t) + J 2 (t) + I 3 (t) + h(t) + I 5 (t) + h(t) + I 7 (t). 

Let us estimate separately the seven terms. From now on all constants depend on 5, 
but are still independent on e. 

First of all we have that I 4 (t) < c 8 G(t) and I 6 (t) = -2F(t). 

Due to the strict parabolicity in the interior and (3.17) we have that 

h(t) < J \ P t\ • K t } 2 dr < c 9 J y?"«)[< t ] 2 dr = c 9 G(t), 

h{t) < cio J W"K)\ ■ \u e „\ • [u e rt ] 2 dr < Cll G(t). 

Moreover from inequality 2ab < ua 2 + u^b 2 we deduce that 

hit) < d 2 J vWRR? • Vv"(K)[Kt] 2 dr < -F(t) + c 13 G(t), 
and similarly 

h{t) < l -F{t)+c u G{t). 

In the same way, exploiting once again the strict parabolicity in the interior, and 
estimate (3.17), we have that 

h < Cis J pV"«)KI • Krt\ dr < -F(t) + c 16 G(t). 
Putting all together we obtain that E'(t) < —F(t) + c 17 G(t), hence 

E{t)+ [ F(s)ds < E(0) + c 17 [ G(s)ds. 
Jo Jo 
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The term E(0) depends on the initial condition only. The integral of G(t) is uniformly 
bounded because of (3.19). It follows that the left-hand side is bounded independently 
on e. Thanks to the properties of p(r, t), this proves (3.20) and (3.22). Finally, from 
(4.15) we have that 

vAOKJ < K t \ + ■ K r \ 2 + + 

From (3.17) and the uniform parabolicity in the interior it follows that 

\u £ rrr \ 2 dr < Cis + Ci 9 J \u £ rt \ 2 dr, 

which proves (3.21). 



4.2 Proof of Theorem 3.3 

Let us briefly sketch the outline of the proof, which is quite similar to the proof of 
Theorem 3.1. First of all we show that the initial boundary value problem has a unique 
solution u £ G C 2,1 (T £ ) H C°°(T £ H Intp(T)). Then we set v =: u £ r) and w := u e rr . Both v 
and w belong to C°(T £ ) fl C°°{T £ fl Intp(T)). It is easy to see that v is the solution of 
the following equation 

v t = -tp"{v)v rr - V "\v)v 2 r - ^-v r + ^ V(r, t) e Int P (T e ), (4.50) 

with Dirichlet boundary conditions 

v(P(to-t),t)=v(-y(to-t),t) = l + e Vte[e,f ], (4.51) 
and initial datum 

v(r,e) = l + e Vr e \J3{t - e), 7 (i - e)]. (4.52) 

In the same way w is a solution in Intp(T^) of equation (once again the terms in v 
are thought as coefficients) 

w t = —Lp"(v)w rr — 3(p"'(v)w r w — Lp IV [v)w z 

_^M w ,_m Wr+2 ^i w . 2 ^ (4 . 53) 

with Dirichlet boundary conditions 

w(P(t -t),t)= u £ rr ((3(t -t),t) Vt G [e, t ], (4.54) 
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w(>y(t -t),t) 



KMto -t),t) 



Vi e Mo], 



(4.55) 



and initial datum 



w(r,e) = Vr e{f3(t -e) n (t -s)]. 

Then we show that v and u> satisfy four sets of inequalities in %. 
• The first pair of inequalities is 

1 + e < v(r,t) < 2 + (p'(l)t. 



(4.57) 



(4.56) 



From the second inequality in (3.4) we have in particular that t < l/<//(l). There- 
fore (4.57) implies (3.26). 

• The second inequality is that 



for a suitable constant rj > 0. The term after rj is positive in Intp(T). This implies 
(3.27), hence also (3.28). 

• The third pair of inequalities is 



v(r, t)>l + e+ (b(t -t)-Vi)(r- /3(t - t)) - 7 o (r - (5{t - t)f , (4.59) 
v(r, t)<l + e+ {b{t -t) + y/i)(r- f3(t - t)) + 7 o (r - f3(t -t)) 2 , (4.60) 



where 70 is the constant defined in (3.1). Arguing as in the proof of Theorem 3.1, 
these inequalities yield (3.29). 

In an analogous way we have that 



v(r, t)>l + e+ {c{t -t) + y/i)(r- 7 (f -t))- lo (r- 7 (i - i)f , (4.61) 
v(r, t)<l+e+ (c(t -t)-y/i){r- 7 (* - *)) + 7 o (r - 7 (*o " t)f , (4.62) 



which imply (3.30). 

• Thanks to (3.29) and (3.30) we have an estimate on the values of w at the bound- 
ary. This is the starting point to prove the fourth set of inequalities 



b(t -t)-^-li(r-(3(t -t)) < w(r,t) < b(t -t) + ^+li(r-(3(t -t)), (4.63) 
c(t -i)->/i-7i(7(fo-*)-r) < w(r,t) < c(t -t)+^+^(t -t)-r) : (4.64) 



where 71 is the constant defined in (3.2). These imply (3.31) and (3.32). 

We are now ready to proceed with the details. Many steps of the proof (for example 
the integral estimates) are analogous to the corresponding steps in the proof of Theo- 
rem 3.1. In these cases we skip them, focussing only on what is different. We also skip 
the proofs of (4.61), (4.62), (4.64), which are analogous to the proofs of (4.59), (4.60), 
and (4.63), respectively. 




(4.58) 
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Existence and maximum principle for space derivatives Let us take a function 
(f £ G C°°(R) which coincides with ip in the interval [1 + e, 3], and such that <p"(o~) < 
—v e < for every a G R. Equation (3.24), with <p £ instead of ip, is strictly forward 
parabolic. Therefore by well know arguments the initial boundary value problem admits 
a unique solution u £ G C 2 '\T £ ) n C°°(T e \ d P (T)). 

It is easy to see that z(r, t) :— 1 + e is a subsolution of problem (4.50) through (4.52), 
both with ip £ and with <p. This proves the lower bound in (3.26). 

We claim that z(r, t) := 2 + <p'(l)t is a supersolution of (4.50) through (4.52). Indeed 
on the three sides of the parabolic boundary of T e we have that 

z(r,t) >2>l+e = v(r,t). 

In the interior we have that z r (r,t) = z rr (r,t) = 0, hence (since 1 < z(r,t) < 3) 

z t {r-t) = V \l)> V \z{r-t))> ^M2 , 

which is exactly the required differential inequality (both with <p e and with <p). 
This completes the proof of (4.57), hence of (3.26). 

Uniform strict parabolicity in the interior Let us choose 7] > small enough so 
that ^ 

V<t (<l), V<Q- Q + 20 ^) ^'(3). (4.65) 

Let z(r,t) denote the right-hand side of (4.58). We claim that, under these restric- 
tions on rj, z is a subsolution of (4.50) through (4.52), which implies (4.58). 

Parabolic boundary When r = [3(t — t) or r = 7(t — t) we have that 

z((3(t -t),t) = z(rf(t -t),t) = l<l + e = v((3(t -t),t) = v(rf(to -t),t). 
Since i] <t , when t = e we have that 

z{r,e) = I + 77 ( - - (^-3) 2 ) < 1 + — < l+e = v{r,e). 

Differential inequality We have that 

z r (r,t) = -2(r - 3)?7, z rr (r,t) = -2rj, z t (r,t) = r)/t Q . 
In order to show that z is a subsolution we have to prove that 

^ < 2 W '\z) - ^\"\z){r - 3) 2 + 2 V ^-(r - 3) + 

to r r 2 
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It is easy to show that 1 < z(r,t) < 3 for every (r,t) G T £ . Therefore from the 
properties of <// we have that (p'(z)/r 2 > (l/16)<//(3). All the other terms are uniformly 
small when 77 is small. This shows that the differential inequality is satisfied when 77 is 
small enough, for example as soon as 77 fulfils the last inequality in (4.65). 

Estimate on second derivatives at the moving boundary — Subsolution Let 

us prove (4.59). To this end we denote its right-hand side by z{r,t), and we show that 
it is a subsolution of (4.50) through (4.52). 

Boundary r = j3(t — t) In this case we have that 

z(/3(t -t),t) = l + e = v(/3(t -t),t). 

Boundary r = j(to — t) From the explicit expressions (2.7) we have that 

z(>y(t -t),t) = l+e+ {b{t - t)-y/i) -2^-470^- 

yto to 

< l + e + 2-^(&(fo-f)-27o^L). (4.66) 

Exploiting estimate (4.13), and the fact that 70 > 1/2 and t < 1, we have that 

\ft 

b(t -t)<Vt< 2 l0 Vi < 2 7o-^=- (4.67) 
From (4.66) and (4.67) we conclude that z(j(to — t),t) < 1 + e = v(j(to — t),t). 

Boundary t = e From (4.13) we have that b(t — e) < y/e, hence 

z(r,e) = l+e+(b{t -e)-^)(r-P{t -e))- l0 {r-P{t -e)f 
< 1 + e = v(r, e). 

Note that the term with ^/e instead of e in the right-hand side of (4.59) is essential 
in this point of the proof. 

Differential inequality Let us set for simplicity x := r — f3(to — t). Then we have 
that 

z r (r, t) = b(t - t) - y/e - 27 x, z rr (r, t) = -2j , 

z t (r, t) = b(to - t)P'(t -t)- P \t - t)y/e - (b'(t -t) + 2 lo p'(t - t))x. 



32 



Plugging these expressions in (4.50), and exploiting the equation defining b(t — t) 
as we did in (4.36), we end up with the following differential inequality 



P'(t -t)y/i+ (b'(t -t) + 2 lo (3'(t -t))x > 
> -2 7o¥ >"(z) 

+<p"'(z) (b(t -t)-y/i- 2 l0 x) 2 - tp'"(l)b 2 (t - t) 

+ < ^(Kt -t)-y/i-*y x) 



(3 2 (t -t) r 2 ' 

This inequality has to be proved for all (r,t) G T £ such that 1 < z(r,t) < 3. This 
can be done arguing as we did in the corresponding step of the proof of Theorem 3.1 
(and using from time to time that -Je > e). 

Estimate on second derivatives at the moving boundary — Supersolution In 

order to prove (4.60) it is enough to show that its right-hand side is a supersolution of 
(4.50) through (4.52). In this case the inequalities on the three pieces of the parabolic 
boundary of % are trivial. The proof of the differential inequality is analogous to the 
case of the supersolution. We skip the details. 

Global estimate on second derivatives — Subsolution Let z(r, t) denote the left- 
hand side of (4.63). We claim that it is a subsolution of (4.53) through (4.56). 

Boundary r — (3 (to — t) From (3.29) we have that 

z(P(t ~t),t) = b(t -t)-y/e< u £ rr (P(t -t),t)= w((5(to -t),t). 

Boundary r = 7(t — t) Exploiting estimates (4.13), (4.14), and (3.30), and the fact 
that 71 > 1 and t < 1, we have that 

z(rf(t - t),t) = b(t -t)-y/E- 2 7l ^l <Vt-Vi- 2 7l ^l < -Vt - V~e < 

yto vh) 

< c(t -t)-y/e< u £ rr (l(to ~t),t) = w(rf(t -t),t). 

Boundary t = e From (4.13) we have that b(to — e) < y/e, hence 

z(r, e) = b(t - e) - ^fe - 71 (r - f3(t - e)) < = w(r, e). 
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Differential inequality After computing the derivatives and changing the signs we 
have to prove that 

6<( t0 - t)+7l/ 3<( t0 - t ) > -3 7l ^>) 2+ ^ Wz = + ^^- 7l ^-2^M,+2^. 

Note that this is exactly the same inequality satisfied by supersolutions in Qi (since 
we reversed the time, subsolutions correspond to supersolutions and vice versa). The 
proof is of course completely analogous. 

Global estimate on second derivatives — Supersolution Let z(r, t) denote the 
right-hand side of (4.63). We claim that z is a supersolution of (4.53) through (4.56). 

Boundary r = /3(to — t) From (3.29) we have that 

z(P(t -t),t) = b(t -t) + ^> u £ r Mh ~t),t)= w(P{t -t),t). 

Boundary r = 7 (i — t) We have that 

Vt 

z( 7 (t -t),t) = b{t - t) + + 2 7l -^= > yfl. 

Moreover, from (3.30) and the fact that c(t — t) < 0, we have that 
iw( 7 (t -t),t) = < r ( 7 (t -*),*)< c (to - t) + < v^- 
It follows that z(7(t„ -t),t)> w( 7 (t — *),*)■ 

Boundary t = e We have trivially that z(r, e) > = w(r, e). 

Differential inequality We have to prove the same inequality satisfied by the corre- 
sponding subsolution in Q 1 . The proof is the same. 

5 Passing to the limit 

The basic tool is the following compactness result. We omit the proof, for which we refer 
to [10, Lemma 3.1]. The key point in that (5.1) through (5.3) yield a uniform bound on 
the norm of f £ in the Holder space C 1//2 ' 1 / 4 ([r 1 , r 2 ] x [ti,t 2 ]). 
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Lemma 5.1 Let [r 1 ,r 2 ] x t 2 ] be a rectangle, and let e > 0. For every e G (0, e ), 
let f £ G C 1 ([r 1 ,r 2 ] x [ti,t 2 ]). Let us assume that there exists Mel such that 




ir(ri,ti)| < M VeG(0,5 ), (5.1) 
[fr(r,t)} 2 dr < M V* G [h, t 2 ], V £ G (0, e ), (5.2) 



[f £ {r,t)} 2 drdt < M Ve G (0,e ). (5.3) 



Then the family {f £ } is relatively compact in C°([ri,r 2 ] x [ti,t 2 ]). 1=1 



5.1 Proof of Theorem 2.1 



Let to an d Uq be as in Theorem 3.1, and let {M e } ee (o 5 i) be the family of solutions produced 
by that theorem. We claim that u e uniformly converges, as e — * + , to a limit u, which 
satisfies the initial condition 



and all the conditions required in Theorem 2.1. 

Uniqueness Let us assume that a function u G C 2,1 (Qi) satisfies equation (2.6) 
for every (r,t) G <2i, the Neumann boundary condition (1.5) for every t G [0,t ], the 
Neumann boundary condition (2.16), and the initial condition (5.4). Let us assume 
also that < u r (r,t) < 1 for every (r, t) G Q,±. Under this condition equation (2.6) is 
degenerate parabolic, hence the solution u is unique. 

This shows that the limit problem provides a unique characterization of the possible 
limits. Therefore in what follows we can limit ourselves to show that u e converges to 
the solution of this problem up to subsequences. 

Convergence of u e Estimates (3.11) and (3.18) provide a uniform bound on the 
Lipschitz constant of u £ . Moreover the functions u £ are equi-bounded due to the initial 
condition. From the classical Ascoli's Theorem it follows that (up to subsequences, which 
we don't relabel) u e uniformly converges in Q\ to a continuous function u, satisfying of 
course the initial condition (5.4). 

Convergence of u\ Let us consider any rectangle [ri,r 2 ] x [ti,t 2 ] Q Qi \ T\, and let 
us apply Lemma 5.1 in this rectangle with f £ = u £ . In the rectangle assumption (5.1) 
follows follows from (3.11), assumption (5.2) follows from (3.17), and assumption (5.3) 
follows from (3.19) and the strict parabolicity in the interior. We obtain that 



u(r,0) = u (r) VrG[l,2] 



(5-4) 



uniformly in [ri,r 2 ] x [ti,t 2 ]. 
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Taking the union over all such rectangles, we easily conclude that 

u e r (r,t) ^ Ur(r,t) V(r,*) e 

We claim that u r can be continuously extended to the whole Q\ by setting u r (r, t) — 1 
for every (r,t) G r\. In order to prove the continuity of this extension, let us consider 
(fi,ti) = (f3(ti),ti) G Ti and (r 2 ,t 2 ) G Qi \ IY From (3.17) we have that 

|<(r,0 - (l + e)| = \u e r (r,t)-u e r (P(t),t)\ < M 4 \r-f3(t)\. 

Letting e — > + we obtain that 

\u r (r,t)-l\ <M 4 \r-P(t)\ V(r,t) G Qi \ I\. 

Therefore we have that 

\ur(r 2 ,t 2 ) - w r (ri,ti)| = K(r 2 ,i 2 ) - 1| 

< M 4 |r 2 -/3(f 2 )| 

< M 4 (|r 2 -ri| + |/3(ti)-/3(t 2 )|). 

Due to the continuity of (3{t), the right-hand side is small when (r 2 ,t 2 ) is close to 
(ri,^i). This proves the continuity of u r up to r\. Of course u r satisfies the Neumann 
boundary conditions in r = 1 and in r = /3(t). 

Convergence of u e rr Let us consider any rectangle [r 1; r 2 ] x [t 1 ,t 2 ] C Q x \r l5 and let 
us apply Lemma 5.1 in this rectangle with f £ = u £ rr . In the rectangle assumption (5.1) 
follows follows from (3.17), assumption (5.2) follows from (3.21), and assumption (5.3) 
follows from (3.22). We obtain that uniformly in [n,r 2 ] X [h,t 2 ]. 

Taking the union over all such rectangles, we easily conclude that 

u e „ (r, t) -> u rr (r, t) V(r, t) G Q 1 \ r\ . 

Let us extend u rr to the whole Q 1 by setting u rr (r, t) = b(t) for every (r,t) G IY In 
order to prove the continuity of this extension, let us consider (ri,ti) = (/3(ti),ti) G Fi 
and (r 2 ,t 2 ) G Qi \ IV Passing (3.16) to the limit we obtain that 

\u rr (r,t)-b(t)\ < M 3 \r-P(t)\ V(r,t) G Q 1 \T 1 , 

hence 

\u rr (r 2 ,t 2 ) - Urr(ri,fi)| = \u rr (r 2 ,t 2 ) - 6(*i)| 

< |Mrr(r 2 ,t 2 ) -6(t 2 )| + |6(* 2 ) -6(*0I 

< M 3 |r 2 -/3(t 2 )| + |6(t 2 ) -6(tx)| 

< M 3 (|r 2 - n| + \0(t 2 ) - Pih)]) + |6(* 2 ) - 6(tx)|. 

The conclusion easily follows from the continuity of and 
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Convergence of u\ Since u\ is related to u\ and u e rr by (2.6), the convergence of 
u\ to a continuous function in Qi follows from the convergence of u\ and u £ rr and the 
continuity of their limits. The limit of u\ is of course u t . This completes the proof that 
u is of class C 2 ' 1 and solves equation (2.6). 

5.2 Proof of Theorem 2.2 

We have already seen in section 3.2 that Theorem 2.2 is equivalent to Theorem 3.2. In 
turn, the proof of Theorem 3.2 follows by passing to the limit the solutions u e provided 
by Theorem 3.3. 

The argument is quite similar to the proof of Theorem 3.1. Thanks to some of the 
estimates of Theorem 3.3 we can apply Lemma 5.1 and pass to the limit in any rectangle 
[r 1; r 2 ] x t 2 ] C Intp(T), where u £ is defined for every e small enough. Thanks to the 
remaining estimates of Theorem 3.3 we can show that the limit u and its derivatives u r , 
u rr , u t can be continuously extended up to the boundary. We skip the details. 
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